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We present and discuss constraints on supersymmetric type I seesaw models imposed by 
neutrino data, charged lepton flavor violation and thermal leptogenesis. 



1 Supersymmetric seesaw mechanism and slepton mass matrix 

The observed neutrino oscillations imply the existence of neutrino masses and lepton flavor 
mixing, and give hints towards physics beyond the Standard Model. For example, the smallness 
of the neutrino masses suggests the realization of the seesaw mechanism involving heavy right- 
handed Majorana neutrinos. The latter violate lepton number and CP, and thus allow for 
leptogenesis. Particularly interesting are supersymmetric scenarios, where the lepton flavor 
violation (LFV) present in the neutrino sector is transmitted to the slepton sector giving also 
rise to measurable LFV processes of charged leptons. 

A minimal model of this kind is obtained if three right-handed neutrino singlet fields vr are 
added to the MSSM particle content. In this model, one can have the following Majorana mass 
and Yukawa interaction terms 1 : 

1 



vfMv c R + v1(Y v L-H 2 , 



(1) 



where M is the Majorana mass matrix, Y u is the matrix of Yukawa couplings, and L and H2 
denote the left-handed lepton and hypercharge +1/2 Higgs doublets, respectively. Electroweak 
symmetry breaking then generates the neutrino Dirac mass matrix ra^, = Yy^H®) where (H®) = 



■usin/3 is the appropriate Higgs v.e.v. with v = 174 GeV and tan/3 



If the mass scale 



Mr of the matrix M is much greater than the electroweak scale, and thus much greater than 
the scale of mp, one naturally obtains three light neutrinos with the mass matrix 



M v = m^M^rriD = YJM~ x Y v {v sin/?) 2 



(2) 



and three heavy neutrinos with the mass matrix = M. In the basis assumed, M is diagonal, 
while M v is to be diagonalized by the unitary MNS matrix U : 



U T M U U = diag(mi, m 2 ,m 3 ), (3) 
U = FcKM(0i2^i3^23^)-diag(e^,e^ 2 ,l), 

6ij being mixing angles, 5 and fa being Dirac and Majorana phases, respectively, and mi being 
the light neutrino mass eigenvalues. 

The heavy neutrino mass eigenstates N, which are too heavy to be observed directly, influ- 
ence the evolution of the MSSM slepton mass matrix: 

m? = ( < m \A +( 5 f{* ) . (4) 

1 V m lR m R /mssm V 8ml R Sm R J N 

It is these flavor off-diagonal virtual effects which lead to charged LFV. Adopting the minimal 
supergravity (mSUGRA) scheme one finds, in leading logarithmic approximation 2 , 

6m 2 L = -^(3m 2 + A 2 )YjLY u , 5m 2 R = 0, 5m 2 LR = -jjr^Aov cos f3Y l Y^LY u , (5) 

where = \n(MQ\j^/Mi)Sij, Mi being the heavy neutrino masses, and mo and Aq are the 
universal scalar mass and trilinear coupling, respectively, at the GUT scale Mqut- 
By inverting (2), the neutrino Yukawa matrix can be written as follows 1 : 

Y v = —^— diag(/Mi, 7M2, v / M3)-«-diag(V^T, v^> V™s)-U ] '. (6) 

V Sill Jj 

Here, a new complex orthogonal matrix R appears which may be parametrized in terms of 3 
complex angles 9i = x« + iyf. 

-C1S3 - S1S2C3 S1S3 - C1S2C3 
R = I c 2 s 3 cic 3 - S1S2S3 -S1C3 - C1S2S3 ) , (7) 

S\C2 CiC 2 

with (q, Si) = (cos 6i, sin 6,1) = (cos x-i cosh yi—i sin X{ sinh yi, sin Xi cosh yi+i cos sinh yi). While 
the light neutrino masses m; and the mixing angles 6ij have been measured or at least con- 
strained, the phases fa and 5, the heavy neutrino masses Mj and the matrix are presently 
unknown. Using the available neutrino data 3 as input in the appropriate renormalization group 
equations, Y v is evolved from the electroweak scale to the GUT scale and then put into the 
renormalization of the slepton mass matrix from Mqut to the electroweak scale. 

2 Charged lepton flavor violation 

The renormalization effects (5) lead to charged LFV either via contributions of virtual sleptons 
in loops such as in radiative decays — ► Ijj, or via real slepton production and decay such as 
in e + e~ — ► — > + 2\\. To lowest order in LFV couplings one has 1,2 

\(5mr) 2 -\ 2 

T(U Ij-y) oc G?m\ ' V _ tan 2 0, (8) 




m 



m characterizing the typical sparticle masses in the loop. Similarly, one finds 



4 



2 12 



-(e + e- - iri+ + 2x?) - ^-^(e+e" - Zr(+ + 2*°). (9) 

m j 1 f 



Consequently, light neutrino data imply interesting constraints on LFV processes through the 
dependence on \(Y„ LY v )ij\ 2 with Y v given by Eq. (6). Conversely, measurements or bounds on 
LFV processes can constrain the fundamental seesaw parameters Mj and R. 

For simplicity, we first consider the case of mass degenerate heavy Majorana neutrinos with 
Mj = Mr. If R is real, i.e. yi = , then it will drop out from the product YjY u in this case 
as do the Majorana phases fa and fa, leaving Mr and 5 as the only unconstrained parameters. 
Fig. la shows the typical rise of Br(ij — ► ^7) with Mj| suggested by Eq. (8) for fixed light 
neutrino masses. Also shown is the impact of the uncertainties in the neutrino data. From 
the present bound 5 Br(/x — ► ey) < 1.2 • 10~ n one can derive an upper limit on Mr of order 
10 GeV. Furthermore, from Eqs. (8) and (9) the uncertainties in the neutrino parameters are 
expected to drop out of the correlation of radiative decays and scattering processes in the same 
LFV channel. This is demonstrated in Fig. lb for the Tfi channel. As can be seen, combined 
measurements of both processes provide decisive tests of the considered scenarios. 
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Figure 1: (a) Br(r — > /try) (upper points) and Br(/i — > 67) (lower) versus Mb in mSUGRA scenario SPSla for 
real R. The light neutrino masses are assumed quasi-degenerate. The mixing angles and the Dirac phase are 
scattered within the full ranges consistent with experiment. The solid (dashed) horizontal lines mark the present 
6 ' 5 (expected future) bounds, Br(r -» fjrf) < 6.8 x 1(T 8 (1CT 8 ) and Br(p -» erf) < 1.2 x 10" 11 (lCT 13 ). (From 7 .) 
(b) cr(e + e~ — > r + fi~ + 2x5) versus Br(r — > fij). The light neutrino parameters are scattered as in (a). The plots 
(from left to right) are calculated in the mSUGRA scenarios C, B', SPSla, G' and P. (From 4 .) 



The above results are rather conservative, since LFV processes will be enhanced if the light 
neutrino masses are hierarchical instead of degenerate, and/or if R is complex. In the latter case, 
LFV observables have rather more freedom since the dependence on the %)i can be as significant 
as the Mr dependence, as Fig. 2 shows. The change in YjY v is approximately 



A R (YjY v ) « Udteg(y/m)(ltfR- l)diag( v ^)C/ t . 



(10) 



Eq. (10) implies two features seen in Fig. 2: (i) Compared to the case of degenerate light neutrino 
masses, the y dependence in the hierarchical case is weaker since the condition 7713 3> mi 2 implies 
that only (R^R — 1)33 = 0{yf) contributes, (ii) Since Eq. (10) is approximately imaginary and 
linear in the yi , non-zero can only increase the observable so that lower limits obtained for real 
R remain unaffected. Even small values of y can enhance a process by orders of magnitude from 
the real R result. Also, in contrast to the real R case, the LFV branching ratios for degenerate 
neutrinos can now be larger than that for hierarchical neutrinos. 



3 Leptogenesis 

In thermal leptogenesis the baryon asymmetry of the universe is generated from out-of-equilibrium 
decays of the heavy Majorana neutrinos (see e.g. 8 ' 9 ). Most models are based on the assump- 
tion of a hierarchical mass spectrum Ml <C M2 <C M3. It is then natural to assume that light 
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Figure 2: Br(k —* Ij'y) versus yi = y for fixed Mr = 10 12 GeV in mSUGRA scenario SPSla for hierarchical (dark 
points) and degenerate (light) light neutrino masses. The parameters are scattered as in Fig. la. In addition, the 

Xi are scattered over their full range < Xi < 2tt. 



neutrino masses are also hierarchical. In this case, the baryon to photon ratio is determined by 
four factors which are briefly explained below: 

Vb = ~«as P hei«/- (11) 

Thy 

The CP asymmetry ei generated in the decays of the lightest of the heavy Majorana neu- 
trinos N± is 10,11 

e _ r(iVi h 2 + z) - r(iVi -^h 2 + I) _ 3 Mi Ej^fMgij (12) 

61 ~ r (Nx -» h 2 + l) + r(JVi -> h 2 + I) ~ 8vr v 2 sin 2 (3 Y,i m i\Ru\ 2 ' 

This relation clearly shows that non-zero imaginary parts of the R matrix elements are necessary 
to generate a CP asymmetry. The efficiency factor kj in (11) takes into account the washout 
of the initial (B — L) asymmetry. A reliable numerical fit for Kf for hierarchical light neutrinos 
in the strong washout regime can be found in 8 . The solar and atmospheric neutrino mass fits 
favor a value Kf = 0(1O~ 2 ). The (B — L) asymmetry is subsequently converted to a baryon 
asymmetry by sphaleron processes. In the case of the MSSM one obtains the conversion factor 
asph = 2% • Finally, one has to take into account the dilution of the asymmetry due to standard 
photon production, described by the dilution factor d ~ ^ in the MSSM. Confronted with the 
observed baryon asymmetry 12 r]B = (6.3 ± 0.3) • 10 -10 , relation (12) implies a lower bound 
on the M\ scale 11 , e.g. if ei > 10~ 6 , then Mi > 4 • 10 9 GeV. Furthermore, to allow for 
thermal production of right-handed neutrinos after inflation, we require Mi < 10 11 GeV, the 
maximum order of magnitude that the reheating temperature can reach without suffering an 
overabundance of gravitinos, whose decay into energetic photons can otherwise spoil big bang 
nucleosynthesis. In the above mass range, the condition to reproduce the experimental baryon 
asymmetry then puts constraints on the parameters x 2 and X3 of the R matrix 13 as illustrated 
in Fig. 3a. With decreasing Mi < 10 11 GeV, the allowed values for x 2 are concentrated more 
and more at sinx2 = 0. A similar behaviour is found for the angle x%. 

As discussed in the previous section for degenerate Majorana neutrinos, experimental limits 
on Br(fi — ► e~f) can be used to constrain the heavy neutrino scale, here represented by the 
heaviest Majorana neutrino mass M3. This is shown in Fig. 3b. The present bound on Br(fi — > 
erf) constrains M3 already to be smaller than 10 13 GeV. Future experiments are expected to 
reach below M 3 = 10 12 GeV. 
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Figure 3: (a) Region in the plane (x2, Mi) consistent with the generation of the baryon asymmetry r]B = (6.3±0.3)- 
1(T 10 via leptogenesis. (b) Br(fi -> e-y) as a function of M 3 in mSUGRA scenario SPSla, for Mi = 10 10 GeV. The 
solid (dashed) line indicates the present (expected future) experimental sensitivity. All other seesaw parameters 
are scattered in their allowed ranges for hierarchical light and heavy neutrinos. 



4 Concluding remarks 

Supersymmetric LFV is very model dependent. As far as the mSUGRA parameters are con- 
cerned, this is further emphasized in Fig. 4, where a scan in the (mo, mi/2) plane is performed 
showing contours of fixed LFV cross section at the ILC and ft — > ej branching ratio. Clearly, 
the strategy to probe LFV considered here will only be applicable once sufficient measurements 
of the SUSY particles' properties have been made. 

In a given scenario LFV observables, particularly in the fie channel, are very sensitive to 
uncertainties in the neutrino parameters. However, correlations of observables in the same LFV 
channel suffer much less from this. As demonstrated in Fig. lb, such correlations can be rather 
strong, and therefore very useful probes. Although less pronounced, relations also exist among 
observables in different LFV channels like rfi and fie. These can be exploited to improve direct 
bounds. Fig. 5 illustrates this for a variety of seesaw parameters, uncovering a rather good 
correlation between the fie and Tfi channels. One sees that the present experimental bound 
on Br(/x — > ej) can be used to improve the direct bound on Br(r — > ^7) by almost 2 orders 
of magnitude. Interestingly, this result does not depend on whether hierarchical or degenerate 
heavy and light neutrinos are assumed. 

For hierarchical Majorana masses, leptogenesis provides additional constraints. For a heavy 
Majorana neutrino mass spectrum obeying 10 10 GeV « Mi <C M2 <C M3 « 10 13 GeV (in 
mSUGRA scenario SPSla), both LFV and leptogenesis are found to be viable. The heaviest 
Majorana mass M3 is constrained from above by the LFV process ft — > e"f, the lightest mass 
Mi from below by leptogenesis. In addition, the orthogonal R matrix encoding the mixing of 
the right-handed neutrinos, which is parametrized by the angles 9{ = X{ + it/i, is constrained 
according to sinx2,3 — and y« < 0(1) as a consequence of the successful leptogenesis condition 
for small Mi and the requirement of perturbative Yukawa couplings, respectively. Finally, the 
remaining parameter x\ can be constrained from the ratio Br(fi — > e^)/Br{r — ► ^7) 13 . Other 
work along this line can be found in the literature (for references see 13 ) and has also been 
presented at this meeting (see e.g. 14 ). 

Acknowledgments 

We thank H. Pas, A. Redelbach and Y. Shimizu for fruitful collaboration. This work was 
supported by the Federal Ministry of Education and Research (BMBF) under contract number 
05HT1WWA2. 



Figure 4: Contours of the polarized cross 
section a(e + e~ — > jj, + e~ + 2\i) at \fs ee — 
800 GeV in the mo-ni!^ plane (solid lines). 
The remaining mSUGRA parameters are: 
A = GeV, tan/3 = 5, sign(/i) = +. The 
beam polarizations are: P e _ = +0.9, P e + 
■ 0.7 For comparison, Br(/x — > e7) is shown 
by dashed lines. The neutrino oscillation pa- 
rameters are fixed at their best fit values, 
the lightest neutrino mass and all complex 
phases are set to zero, and the degenerate 
Majorana mass scale is Mr — 1 14 GeV. The 
shaded (red) areas are forbidden by mass 
bounds from various experimental sparticle 
searches. 
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Figure 5: Br(r — > ^7) versus Br(^i — > &y) in 
mSUGRA scenario SPSfa, assuming degenerate 
Majorana masses and hierarchical (diamonds) and 
degenerate (stars) light neutrino masses, with pa- 
rameters scattered as in Fig. 1. In the case of hierar- 
chical Majorana and light neutrino masses (squares) 
all seesaw parameters including the Xi are scattered 
within their experimentally allowed ranges, while 
the yi and Mi are scattered within the bounds de- 
manded by leptogenesis and perturbativity. Also in- 
dicated are the present experimental bounds 5,6 . 
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